00 

in 



Consistent BCS+RQRPA formalism with application to the 

double beta decay 

A. Bobyk, Wieslaw A. Kaminski, P. Zar§ba 
Department of Theoretical Physics, Maria Curie-Sklodowska University, 
Radziszewskiego 10, 20-031 Lublin, Poland 
(February 9, 2008) 

A new consistent analysis of the renormalized proton-neutron quasiparti- 
cle random phase approximation based on the simultaneous recalculation of 
the one-body density matrix and the pairing tensor has been used to study 
the double beta decay. We demonstrated that inclusion of the quasiparti- 



> . 

■ cle correlations at the BCS level reduces the ground state correlations in the 
• particle-particle channel of the proton-neutron interaction. We also simpli- 

o ■ 

OO ■ fied the RQRPA equations significantly obtaining a low-dimensioned set of 

linear equations for the quasiparticle densities. The formalism was applied to 

■ the double beta decay in ^^Ge. 

> 

^ ! The proton-neutron quasiparticle random phase approximation (pn-QRPA) has been 

H ; 

5^ , considered the most powerful method for the beta and double beta transition calculations of 
nuclear sytems which are far away from the closed shells [|1],||,0,0,|^J|,0,|| . A remarkable suc- 
cess was achieved especially by the QRPA approach in revealing the suppression mechanism 
of the neutrino accompanied double beta decay, a long-standing problem of the theoretical 
treatment of this process. Further development of the approach went beyond many short- 
comings and refined calculations of nuclear matrix elements involved in the double beta 
decay. Among others, the following problems were set and solved: particle number non- 
conservation P,pi]|, role of the proton-neutron pairing [|rT|, violation of the Pauli exclusion 



principle ||T^,|T3[ , higher-order corrections to the ordinary QRPA [|T^,|r^ , treatment of tran- 
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sitions to final excited states [p^|T7tp^pO[ , extension of tlie definition of plionon operators 



by means of the-so-called scattering terms [21|, etc. 



Most of such improvements disregarded, however, the main source of the formalism in- 
stability connected with violation of the Pauli exclusion principle by using the commutation 
relations for the QRPA phonon operators. To overcome this shortcoming of the pn-QRPA 
framework the renormalization technique was proposed and extended to include the 
proton-neutron pairing |jl3|. The main goal of the method called in literature the renormal- 



ized QRPA (RQRPA) is to use a self-iteration of the QRPA equation to take into account 
the additional one-quasiparticle scattering terms in the commutation relations. But this 
procedure results in a non-vanishing quasiparticle content of the ground states and followed 
by some inconsistency between RQRPA and the BCS approach since the ground state ap- 
proximated by the BSC state is chosen to be the quasipatricle vacuum. To minimize the 
influence of such a discrepancy one needs to reformulate the BCS equations in a way pro- 
posed in . Combining both RQRPA and so modified BCS one obtains the self-consistent 
BCS+RQRPA approach (SRQRPA) which we study in more detail in this paper. 

In the QRPA (either ordinary or renormalized) approach one assumes the harmonicity 
of the nuclear motion and starts with the excited-state creation phonon operators of the 



form 22,24 



Qj^M — X! [-^(^n)J^^lpn)J^M ^(^n)J^^{pn)J''M , (l) 
pn 

where X'^^^^ and y(^n)j-^ forward-going and backward-going variational amplitudes, 

respectively. The quantities ^(p^yTr^^ = ['^pQ-nlj'^M are the angular-momentum coupled two- 
quasiparticle creation operators. Since they do not fulfil the bosonic commutation relations 
exactly in the quasiboson approximation (QBA), that is used to derive the usual QRPA 
equations, the Pauli principle is violated. To avoid this serious drawback in the improved 



version of the theory one introduces the renormalized operators ||19|| : 
along with the renormalized amplitudes: 



vm — r)l/2 V'm 



(3) 
(4) 



where Dpn matrix is defined by the expectation value in the RPA ground-state of the 
commutator: 



0) = (1 - rip - fin) 



(5) 



together with the quasiparticle densities {ja = V2ja + 1): 

'^P = ip^(o|K«p]oo|o), (6) 
^n = r'(o|[ata„]oo|o). (7) 

The above equation has been derived using the exact fermionic commutation relations and 
thus it goes beyond the ordinary quasiboson approximation. One can prove easily that now 
the following relation holds: 







>MM' 



(8) 



i.e. the renormalized operators behave as bosons, at least in the sense of the ground-state 
expectation value of their commutator. The phonon operator now reads: 



Qlhl — X! ['^(pn)J^-^lpn 

pn 



and using e.g. the equation of motion (EOM) method one gets the RQRPA equations 
of the usual form: 



■\-)m 



Qva 



-\}m 



(9) 



with the new renormalized RPA matrices A and B defined in [|T2| . Here Vl^-n = E^n — Eq is 
the RPA excitation energy with respect to the ground-state. 

Now there appears the question of the calculation of the Dpn matrix entering the expres- 



sions for the RPA matrices A and B. Using the mapping ||12|| : 



J^Mn 



(10) 



J^Mp 

and inverting ([l|) one derives the following equations for the quasiparticle densities: 



(11) 



fir, 



j;'J' E Dpn\y, 

J^mn 

Jn'j' E Dpn\y[ 
J'^mp 



m 1 2 

{pn),J'^ I 5 



m 1 2 
pn)J^ I • 



(12) 
(13) 



Insterting into (|T3D one gets the following system of linear equations for Up and n„: 



y'p'^p + E ypni^n = y^ 

n 

y'n^n + E ypn^p = 3^, 



pj 



(14) 



or, in the matrix form: 



where 



diag(3^.' 



ip)) 



y 



ciiag(3^;„ 



( \ 

^ip) 



y(.p) 
^y{n) J 



[W) 



ypn — E I3^(pn)j'^p5 



yp^Y^ypn^ y'p^fp + yp. 

n 

y,^Y.ypn, y'n^fn + yn. 



(15) 



It is worth mentioning, that the dimension of our linear problem is only 2n x 2n, where n 
is the dimension of the single-particle basis. This is of much advantage, since (|14[) has to be 
solved many times as one should iterate between (^ and (^^ until convergence is achieved. 
The other way round, inserting (|13]) into (^, as it has been done by several authors, e.g. 



12,13 , one obtains the equation: 



Dpn 1 jp ^ ] Dpn'ypn' Jn ^ ] Dpinyp'n^ 



(16) 
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that, on the contrary to the claim expressed in |]T2|, can be transformed into the x 
hnear system: 

Wpn,p'n'Dp'n' = U, (17) 

p'n' 

where 

^^pn,p'n' — Spp'^nn' ~l~ ^pp'Jp ^pn' ~l~ ^nn'Jn ^yp'n, (18) 

and U is the vector of I's. In practice however, it takes much less time to solve this system 



using the standard linear algebra procedures than following iteration methods, as in |T2|,|13 
But further reduction of the complexity of the problem to the form of (|1^ allows us to 
take all possible multipolarities into account for the calculation of the Dpn renormalization 
factors. Although some of them are less important and neglected in |T2| , p!3[| , when one 
leaves only a few, the J-coupling scheme breaks, since the basis becomes incomplete and 
the validity of the mapping (^)-(|TT|) is questionable. As can be seen further the results are 
the evidence for it. 

With a non-vanishing quasiparticle content of the ground state one arrives at the incon- 
sistency between RQRPA and BCS, since in the latter assumes the ground state to be the 
quasiparticle vacuum. One thus needs to reformulate the BCS equations p3|, namely by 



recalculating the density matrix p and the pairing tensor k. With the standard Bogoliubov- 
Valatin transformation p4| they read now: 



Pa = (O \cic^\ O) =vl + {ul - vl)na, (19) 

Ka = {0\CaCa\0) = UaVa{l - 2na), (20) 

and depend on the quasiparticle densities n^, where a runs over proton or neutron indices. 
The u and v coefficients are obtained by minimizing the ground-state energy, that by virtue 
of the Wick's theorem is expressed as p4|,^ : 







H 



O) = T^fa^aPa + \j2 Jajb {iaafj^lo+\V\ibbfj7l,+ ) K^K, (21) 
a ^ ah 



with the particle-number constraint: 

No = {o N 0) = ^iX- (22) 



In the above, are the single-particle energies and <^(aa)Jrio+l^l(^^)j^io+/ matrix 
elements of the two-body interaction. 

To solve the SRQRPA equations we start with the ordinary BCS equations, putting 
Up = Tin = 0, than to proceed with the corresponding RQRPA problem (inner iteration), that 
gives us new quasiparticle densities and loop with them back to BCS until the convergence is 
achieved (outer iteration). We arrived thus at the doubly-iterative problem and the question 
of efficient and accurate getting through all the calculation steps becomes very impotrant. 
We then stress again, that without showing that the problem of calculating the D-matrix 
can be reduced to the linear system (|T^) of acceptable size the realization of this task would 
be hardly possible. 

The ground-state to ground-state 2z//3/3 Gamow-Teller matrix elements are expressed 
as follows: 



1 r<--.rr, , r^m' . ^ / A ry . ^ \ ^ I A r7\\ ' ^ ' 



(0+g3(A^ + 2)||tTr+| 


lm')(lm' 






|otgs(A^)) 


1 

2 


n^+ + fi-; + g^-(A^ + i)- 







where the charge-changing transition densities are: 



(o^gsik^+iiC') = e(pii^ii^) (24) 



pn 



(lmlk-r+l|Oi!gs) = E(^'ll^ll^)(^P^"A'(^n)l+ + ^P^n3^(pn)l+)V^pn (25) 
pn 

and the overlap of intermediate excited states is assumed to be expressed as: 

(lm'|lm) = E ('^(p'n)l+'^(™«)l+ ~ 3^(^n)l+3^(™n)l+) • (26) 



pn 



In the above the non-primed (primed) quantities result from the SRQRPA calculations 
based on the initial (final) ground-state. 

To illustrate the differences between the QRPA, the RQRPA and the SRQRPA and to 
show much better stability of the SRQRPA solutions, the results of our calculations as a 



function of the particle-particle (gpp) and particle-hole (gph) factors, renormalizing the bare 
two-body NN interaction are plotted in Figs. 1-3. This commonly used renormalization 
is necessary due to tne nucleus finite size (the bare NN matrix elements are calculated for 
the infinite nuclear matter) and due to the limited dimension of the single-particle basis. 
In our calaculations we used the two-body matrix elements calculated from the Bonn-B 
nucleon-nucleon one boson exchange potential. The single-particle energies are calculated 
from the Coulomb-corrected Woods-Saxon potential with the Bertsch parametrization. We 
used several sets of single-particle levels to see how the choice of the basis influences the 
results. We find weak dependence of the RQRPA and the SQRPA results on the dimension 
of the single-particle basis. On the other hand, the QRPA shows no stability on the chosen 
basis. The conclusion is that the most suitable basis for the calculations of 2zy/3/? decay 
in ''^Ge consists of 16 levels with ^^O as a core. Therefore we used this basis in all the 
further studies described below. To compare with the experiment we have adopted the 
experimental half-life of = (1-42 ± 0.03 ± 0.13) x 10^^ yr from the latest measurement 
by the Heidelberg-Moscow jSP cooperation [^. 

In Fig. 1 the calculated double Gamow-Teller matrix elements Mq^^ as a function of gpp 
for two different (7ph values and for three different QRPA approaches are plotted. We would 
like to stress that in these calculations in the self-consistent iterations of the RQRPA and 
the SRQRPA all the intermediate multipolarities were present. The comparison between 
the QRPA, the RQRPA and the SRQRPA results in the physically acceptable region of 
the gpp parameter 0.8 < gpp < 1.2 shows two main features of the renormalized QRPA. 
First, the inclusion of the ground-state correlations beyond QRPA does not only improve 
the agreement between theoretical calculations and experimental data but also causes the 
stabilization of the dependence of as a function of gpp. Second, the iteration procedure 
for quasiparticle densities which causes treating of RQRPA and BCS on the same footing 
stabilizes the results even further. 

Some authors claim that only the limited set of these multipolarities plays a role in the 
evaluation of the double Gamow-Teller matrix elements |1T2| , [T3| . In Fig. 2 there are shown 
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the results of calculations of Mq^^ as a function of Qpp for a different number of multipolarities 
for the RQRPA and the SRQRPA, respectively. The basis is the same as in Fig. 1, but the 
value of giph parameter is fixed to 1.0. It can be seen why the inclusion of all multipolarities is 
essential to obtain the rehable predictions of the RQRPA and the SRQRPA calculations. The 
solid line in Fig. 2 represents the QRPA calculations, the dot-dashed line the calculations 
with only 1+ multipolarity, the dotted line with multipolarities up to 2+, the dashed line up 
to 3~, the long-dashed up to 5~ and the thick solid line all considered multipolarities up to 
11+. The inclusion of higher multipolarities causes the shift of the collapse of the RQRPA 
and the SRQRPA beyond the value of gpp = 1.0. The additional advantage of the SRQRPA 
solutions is that the calculated matrix elements are less dependent on the Qpp parameter. 

In Fig. 3 the effect of including more multipolarities J'^ in the RQRPA and the SQRPA 
calculations of fixed gp^ = 1.0 is shown. The filled symbols represent the calculations 

for gpp — 0.8 and open symbols for gpp — 1.0. One can see the saturation effect for higher 
multipolarities. The explanation of this behaviour is that the higher multipolarities are 
less collective. Their contribution to the ground-state correlations is much smaller than 
the lower ones. An interesting feature is the virtual independence of the matrix element 
on the number of multipolarities in SRQRPA around gpp = 1.0. It is reflected in Fig. 2, 
where one can see that RQRPA "diverges" when going with gpp from to 1, while SRQRPA 
"converges". In our opinion, this is a clear evidence that the self-consistency between BCS 
and RPA, i.e. between the ground-state and excited-state properties is being restored in 
SRQRPA. 

Finally, we would like to address the question of the Ikeda sum-rule violation. It is well 
known, that in the usual QRPA the Ikeda sum rule is conserved exactly if all the spin-orbit 
partners of the single-particle orbitals are present in the basis, i.e. 

5--'5+ = E|(0gslk^+llOr-E|(0gJk^-l|l^)r = 3(iV-^). (27) 

m m 

The violation is marginal even if one or two of these partners are omitted. But this is not 
the case in RQRPA, where the Ikeda sum rule is violated up to 20% and is similar, up to 



25% in the SRQRPA. There is some difference between two nuclei under consideration, i.e. 
germanium, where the situation gets worse when the ground-state correlations (SRQRPA) 
are taken into account, and selenium, where the results are slightly improved. One can 
conclude, that there is more to this than the ground-state correlations to restore the Ikeda 



sum rule. As already pointed out PB|, the scattering terms present in the /3-decay operators 



can be responsible for this effect |^9|]. They give no contribution to QRPA, because they 



are of the quasiparticle-quasihole character, but they should be taken into account when 
the ground-state quasiparticle densities are not assumed to be zero, like in the RQRPA or 
SRQRPA. It is necessary to extend the form of the QRPA phonon operator (|l]) by including 
new excitation modes, the so-called 5-modes ||21|| . 

In conclusion, we have developed a new method of the nuclear matrix-element calcu- 
lations for the neutrino accompanied double beta decay to the ground state in a frame of 
the consistent BCS+RQRPA approach. Using the 2z^/5/3-decay in germanium ^^Ge as an 
example we demostrated that the inclusion of the ground-state correlations beyond QRPA 
causes the stabilization of the dependence of the Gamow-Teller nuclear matrix element, 
but also weakens their influence on their magnitude because of the additional change of the 
quasiparticle densities during the iteration procedure with the modified BCS solution. 

Unlike ortodox QRPA which requires fine tuning of the Qpp parameter describing the 
particle-particle interaction-strength as the matrix element collapses near the physical 
strength, RQRPA and BCS+RQRPA give stable matrix elements over the whole range 
of physical strength and thus the latter approaches allow for more predictive power than the 
old method. 

Owing to the development of the method to calculate the Dpn normalization factors we 
could take all possible multipolarities into account. Then we were able to avoid possibile 
J-scheme breaking because, by neglecting some of the multipolarities, the basis becomes 
incomplete and the vahdity of mappings (p^)-([Tl|) is questionable. 
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FIG. 1. Dependence of the double Gamow-Teller matrix element Mq^ on the renormalization 
factor in the particle-particle channel gpp. Two different cases for gpp = 0.8 and 1.0 and for three 
different QRPA approaches are shown. Magnitude of the experimental estimate is marked by two 
parallel dotted hues p6| . 



FIG. 2. Influence of the different number of the added multipolarities on Mq^ for two types of 
the renormalization: RQRPA and SRQRPA. Full range of the strength factor gpp is scanned for 
fixed gph = 1.0. SRQRPA shows less dependence on gpp parameter than the RQRPA approach. 



The experimental values |26| are between two parallel dotted lines. 



FIG. 3. Contributions of the individual multipolarities (1^, 0^, ... , 11^) to the double 
Gamow-Teller matrix element within two different RQRPA approaches. Calculations were made 
for gph =1.0 and for two different values of the renormalization constant in the particle-particle 
channel: gpp = 0.8, 1.0. For details see the text. 
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